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SUMMARY 



Computations have been made to determine .the effect 
of changes in wing aspect ratio, additional side area, 
flight-path angle, and normal acceleration on the relation 
"between the fin area and the dihedral angle required for 
spiral and for oscillatory lateral stability for a hypo-* 
thetical airplane of the pursuit or fighter categories. 
The calculations, however, are applicable to any type of 
airplane characterized by the parameters and the data em- 
ployed. The results of the computations are presented in 
the form of diagrams of fin area plotted against dihedral 
angle showing the combinations of these variables for which 
the airplane has both spiral and oscillatory stability. 

The diagrams indicate that the effect of wing aspect 
ratio on lateral stability is small. The increased cross** 
wind force that accompanies an increase in side area is 
advantageous in that it makes less difficult the selection 
of a fin area and a dihedral angle that will give lateral 
stability at high speed as well as at low speed with flap 
extended, The effects of flight path and normal acceleration 
are such that airplanes, which are laterally stable in level 
flight, may become unstable in a glide or a climb or when 
subjected to normal acceleration* 



INTRODUCTION 



The results of investigations reported in references 1 
and 3 show that the present-day trends in airplane design 
toward higher wing loadings and larger values of radii of 
gyration in roll and yaw make the attainment of lateral 
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stability difficult , if not' impossible. It has become 
important, therefore, to study in detail all factors that 
affect the lateral stability of the airplane. The limits 
of fin area and dihedral angle required for lateral stabil- 
ity of an airplane as affected by variations in airplane 
density, radii of gyration in roll and yaw, wing chord, 
and tail length have previously been computed and presented 
in reference 1, 

The present report is a continuation of the work of 
reference 1 and includes the effects on lateral stability 
of changes in aspect ratio , side area, flight-path angle, 
and normal acceleration. The results are given in the form 
of diagrams of dihedral angle against fin area for neutral 
spiral and oscillatory divergence* 



A I BP LAKE PARAMETERS 



The airplane considered in the investigation was 
assumed to have the following characteristics: wing area 
S w , 200 square feet; total weight Yf r 6000 pounds; wing 
loading W/S, 30 pounds per square foot; ratio of radius 
of gyration in roll to wing span k^/b, 0.125; and ratio 
of radius of gyration in yaw to wing span kg/b* 0.175. 

The changes in the parameters studied are given in 
the folio wing table: 
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where 



ratio of airplane desnity to air density computed 
for standard air density at sea level (m/pS w b) 

m airplane mass 

p. air density 

S w wing area 

b wing span 

A aspect ratio 

Bf fin area 

Sf additional side ar§a 

a 

l^/b ratio of distance from center of gravity of airplane 
to center of pressure of fin to b 

l a /b ratio of distance of ASf ahead of center of 
gravity to ' "b a 

y flight-path angle 

g normal acceleration 

For the computations, one parameter was varied at a 
time while the others were kept at the mean value. These 
mean values were the same values that were used for like 
quantities in reference 1* 

With regard to the variation of aspect ratio, the 
changes in ^ and t t /b result from the change in b 
that accompanies the change in Ao The wing loading and 
the value of \^ are constant for all values of A. Prom 
elementary considerations the ratio of changes in radii of 
gyration to wing span in roll kx/h and in yaw k z /h would 

be expect ed to change with A J an examinat ion of exist in g 
data from present-day types of airplane, however, showed no 
general tendency for these values to vary with A# 

The additional side area Sf a is used to represent an 

increase in the cross-wind, force of an airplane over that of 
the streamline fuselage for which the original data were 
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obtained. For convenience, the additional side area Sf 
is divided into two parts: ASf a t which is added ahead of 
the center of gravity, and A Sf , which is added to Sf • 
The value of ASf is so chosen that it will exactly counter- 
act the unstable yawing moment due to sideslip of ASf a ; 
that is, the weathercock stability of the airplane is not 

affected by changes in Sf • 

a 

Because few data are available for some of the aero- 
dynamic characteristics required, the lateral-stability 
derivatives were computed from the dimensions and parameters 
given, The contributions of each component of the airplane 
are added to give the values for the complete airplane* 

In the computation of the lateral-stability deriva- 
tives, the changes in the aerodynamic characteristics due 
to changes in the interference effects between the component 
part of the airplane and the effect of changes in power are 
assumed to be zero 0 Actually, for an airplane, some of these 
effects may be large, especially for conditions other than 
straight level flight. Some of the interference effects 
caused by wing location on the variation of the derivatives 
that depend upon sideslip are given in references 3 and 4, 
It follows that the derivatives used in the calculations 
may not represent the actual characteristics of the geometric 
arrangement assumed. The results, therefore, should be re- 
garded as applicable to an airplane characterized by the 
stability derivatives used rather than by a given geometric 
arrangement » 

The stability derivatives used and their variations 
with the airplane parameters were computed by the following 
relations, which are intended to include the effects of the 
parameters and particularly the changes in fin area, addi- 
tional side area, and dihedral angle* These relations are 
the same as given in reference l r modified for constant wing 
area and with a term added to include additional side area. 
The stability derivatives -are the instantaneous rates of 
change 'of the aerodynamic coefficients with attitude or 
angular velocity when the attitude angle or the angular ve- 
locity is zero. For convenience, the partial derivatives 
are written in the form 



dC, ac 1 ac l 

C, Ci (J, etc. instead of * , etc. 

2V 21 
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where 

C ^ force or moment coefficient (rolling-moment coeffi- 
cient in this case) 

P angle of s ides lip 

p rolling velocity 

r yawing velocity 

V velocity of f 1 ight 

The relation for the derivatives is 

Cy = Cy + Cy + Cv + Cv 

x 0 ^(fuselage) X P(S W ) X 3(S f +AS f ) X p(ASf a ) 

(Sf+ASf) ASf 
« -0.16 + 0 - 3.48 ~ — — - 3.48 — -ift (1) 

s w s w 

where 

Cy force coefficient along Y axis of airplane 

p angle of sideslip, positive when the' right wing is 

into the wind, radians 



and the subscripts indicate the contribution of the corre- 
sponding part of the airplane to C v ♦ 

*P 

The constant 0,16 in Cy was computed from 

P(fuselage) 

the relation given in reference 1 for an aspect ratio of 8 
and represents the condition for a constant ratio of wing 
area tc fuselage size, 

Although Cyp(s w ) varies with dihedral (references 3, 

5, 6, and 7), this derivative has not been included because 
it is counteracted, at least partly, by the derivative of the 
side force due to rolling Cy • 
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The rate of change of normal force on the fin and 
additional side* area with angle of sideslip, Cyp ( Sf+ASf ) an<i 

(Sf+AS-p) ASf 
°Y§(Sf )* is ^^a 1 t0 -3.48 — i 5-L and -3.48 & 

a S W S W 

in terms of the wing area* The term Sf/S w is the ratio 
of fin area to wing area (a variable for this investigation) 
and 3o48 is the rate of change of normal-force coefficient 
on the fin with sideslip angle [3 for an effective aspect 
ratio of 3* The value obtained with the fin used on the 
model in reference 3 is about 3*48, 



rolling- moment coefficient 
Z X T 

varies with aspect ratio; it is equal to -0.0175 
for A = 16, -0.0141 for A - 8, and -0.0114 
for A = 4. The values of Ki were obtained 
from data given in reference 6 for a wing of a 
2:1 taper ratio. 

effective dihedral angle of wing, degrees 

Effective dihedral angle is used throughout this report 
as a fictitious angle that would give the wing the value of 
C i p • Wing plan form and elevation, as well as large inter- 
ference effect s , contribute various amounts of Oig. (See 
references 3, 5, and 7.) 

The contribution of the fin to the value of Cin is 



where 
K-i 



r 



fx It Y (Sf+ASf A 

VJ ~ T sini^3.48 — ~ J 



where sin i] is the assumed ratio to the wing 

\b b / 

span of the vertical distance of the center of pressure on 
the fin above the center of gravity* The value of 'x was 
assumed to-be the distance from the. fuselage center line to 
the center of pressure of the fin* Upon the assumption that 
the distance from the top of the fuselage to the. center of 
pressure of the fin varied with, the square root of the fin 
area and by the use of data from reference 3, the relation 



^ » 0#025 + 0,23 V (Sf+ASf ) /S w was obtained and used for 

all variations of the fin, The value of i t the angle of 
fuselage center line to the horizontal, was zero for all 
cases when the lift coefficient Cj, was equal to Q,2« ?or 
other values of C^ f the value of i depended upon "bpth 
C L and A, The value of i for Oj, " 3.8 was the same 
as for Cj, » 1,4, The increase in . was obtained by 

merely deflecting the flaps with no change in the angle of 
attack. The term 3,48 (Sf+ASf)/S w was used in the express 
si on for Cy n , * M - \ 

The contribution of the side area added ahead of the 
center of gravity is 



which is similar to that of Sf + ASf except that tho cen- 
ter of pressure of ASf was assumed to be on the fuselage 
oenter line. a 



-0,204 l t 
— - * K a + rr- 



where C n is the y awing-moment coefficient, 

A value of K a of 0,009 was used for plain wings and of 
0,030 when the flaps were deflected. Values for K 3 may 
vary considerably with wing forms, (See references 5 and 7,) 



The value C n p ( fusela/re ) = -0.204//^ was obtained 

from data in- reference 3.- The change from the expression 
given in reference 1 -was necessary- "because for this in- 
vestigation the fuselage and the wing area are constant. 

The contribution of the fin and the side area to the 
value of 0 no is the computed variation of the normal 

P 

force with p on the fins c Y p(Sf+ASf ) and ° T 3.(AS f ) ; 

times the nondiiaensi onal lever arms l t /l) and l a /T) f re- 
spectively. 

% ' ° l P<s„) " K * U) 

where the values of K 3 (from reference 60 are -0,6 for 
A = 16, -0.5 for A = 8, and -0.4 for A = 4. The 
values of • C lp(Sf+ASf) and * l p ( f uge lage ) are probably 

very small as compared with 0 l p ( w / ag ) and therefore have 
not "been used. 

° n P = %(s w ) + ° n p(s f +As f ) + ° n pUSf a ) 

= k 4 c l + - - 4r sln 0 I 3 - 48 - % ; 

* o V b o J \ S w / 

. siiCla sin i") f 3.48 (5) 

where K 4 varies with the aspect ratio and is (from ref- 
erence 6) -0.089 for A = 16, -0.065 for A = 8, and 
-0.040 for A = 4. 

The expression for °n p ( g f + As f ) ' that is » 

2—2-1 — - — - sin i j i o.48 t_ 

b \b b / ;. S w J 

is the rate of change of the yawing moment with the rolling 
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velocity, produced "by assuming that the normal force on 
the fin is proportional to the angle induced at the cen- 
ter ,pf pre s sure, pf .the fin Jby.,_t.h© rolling velocity; ac- 
tually, the induced angle is a variable along the fin # 

The expression for °np(^s f ) is the same as , for 

a 

Gn p(Sf + A5f) the center of As f a on the fuselage 

center line . 

° l r = ° l r(S w ) + ° J r(S f+ AS f ) + °* r (A8 f ) 

a 

l t fx U \ f S f +AS f \ 

- 2 ia(^ .i. i)(s.«e^) (6) 

where the value of 0.25 was obtained from data given in 

reference 6, and the terms for Ci . . and Ci /A _ x 

^r(S f +AS f ) *r(ASf ) 

a 

are of the same form as for C np(Sf+ASf) and ^^ USf )i 

a 

"but, in this case, the angle is induced "by the yawing ve- 
locity and the values are probably more nearly repre- 
sentative of the actual values of the derivative than the 
values given for C np(Sf+ASf) and ^r^^y 

° n r = ° n r( fuselage) + ° n r(S w ) + ° n r(S f +AS f ) + ° n r(AS f ) 



-0.0242 

+ 



Si 

■ a - AS 



-<v) ( B -«-e) (7) 



where K 5 (from reference 6) is -0,0113 for A « 16, 

-0.02105 for A « 8, and ,-0. 03838 for A = 4 . The con- 
stant K 6 depends on the profile-drag coefficient of the 
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wing (reference 6) and is assumed to be 0,003 for the plain 
wing and 0.030 for the wing with flaps deflected..' The value 
-0,0242/yX was computed from the relations, given in ref- 
erence 1 for the fuselage. This form is used' in th.e- pres- 
ent investigation "because the ratio of wing area to fuse- 
lage size is constant and °n r ( f use i age ) varies inversely 
with span . 

The expressions for °n r ( Sf+ASf ) and c n r(ASf j - 

that is, »(2t) ( 3 .4e^ S -£) and 2 (^) (s.«8^) - 

are the nondimensi onal forms of the rate of change of yaw- 
ing moment due to the fin with yawing velocity. 

CALCULATIONS 

The "boundaries of neutral spiral and oscillatory sta- 
bility wore computed "by use of the lateral- stability equa- 
tions from the theory of small oscillations as given in 
reference 8* Lateral stability depends upon the values 
and the algebraic sign of the term E of the stability 
equation and South 1 s discriminant. When the value of 3 
becomes negative, the airplane becomes spirally unstable. 
The lator&l oscillations increase in amplitude when the 
value of kouth : s discriminant becomes negative. ■ The lim- 
its of the stable region are therefore defined by the 
values of r and Sf/S w when B = 0 and Routh's discrim- 
inant = 0. 

RESULTS AO DISCUSSION 

The results of this investigation are presented in the 
form of diagrams, figures 1 to 4, showing the variations 
in the computed boundaries of spiral and oscillatory sta- 
bility with Sj/o w , the ratio of .fin area to wing area, 
and with r, the effective dihedral angle. The values 
of Sf/S w where C n ^ is zero and 0.05 and the values of 

C t when F is 10° are indicated on the figures. The 

weathercock stability is neutral where C n is zero.- 
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The results, in general t indicate that the value of 
T required for spiral stability increases with 0 L and 
s f/ s w T3aQ value of S f /S w required for oscillatory sta- 
bility decreases with r for several degrecs : and then 
incroasos with continued increase in T. Tor some cases 
the diagrams do not. cover a sufficient range of T to 
show the increase in Sf/S w . TJie rate of change of Sf/S w 
with, r required for oscillatory stability increases with 
C L- ; ■ 

The' stability boundaries for Cj, = 2.8 may appear to 
be inconsistent with the other lift coefficients. The 
differences, however, result from the changes in Cj J9 
C n ^ , and produced by the flaps and from the assump- 

tion that the- angle of attack is the same as for the con- 
dition with no flaps, that is, C L = 1.4. 

Effect of changes in aspect ratio on stability bounds 
a r i e s . - lor spiral stability,' increasing A increases 
the minimum permissible value of F and decreases the al- 
lowable value of S f /S w . (S-ee fig. 1.) For oscillatory 

stability, increasing A decreases the value of Sf/S w 
required for the normal range of values of r, but for 
very large values of F at ' the high lift coefficients, 
this effect may be reversed. The resulting effect of in- 
creasing aspect ratio- on lateral stability is to reduce 
slightly the required value of S f /S w and to increase the 
required value of T. 

E ffect of changes in side area on stability bound- 
arie s*- Adding side area as assumed in this invest!- 

A a 

gation increases the damping coefficients . Cy^ and C Ur 

nnd has no effect on the weather cock-stability coefficient 

^n fl • Increasing Sf /S w decreases the minimum permissi- 
p a 

ble value of F for spiral stability and allows larger 
values of S f /S w to be used. (See fig. 2.) ,Por oscilla- 
tory stability increasing S f /S w ' reduces the required 
value of Sf /S w , 

The resulting effect of increasing . Sf ./S w .. is- to al- 

a 

low much greater variations in the value of Sf/S w that 
can be used. Another variation of side area has been re- 
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ported in reference 1, that is, the case of changes in 
tail length. Increasing the tail length through the nor- 
mal range and decreasing S f /S w so that 0 n ^ remains 

constant produces only minor changes in the stability 
"boundaries. Decreasing Sf/S w and increasing the tail 
length increases C n ^ "but decreases the side area and 

therefore Cy^ . Adding the side, area increases both Cy^ 

and C n . It appears that, in critical cases, it would 
r 

"be "best to add side area as far from the center of grav- 
ity as possible* 

Effect of changes in flight -path on the stability 
boundaries A change in the flight path from a dive to a 
climb greatly reduces the allowable range of Sf/S w re- 
quired for spiral stability. (See fig* 3.) 

For the range of values of T normally used, flight 
path has very little effect on the oscillatory stability 
boundaries. At the high lift coefficients and large val- 
ues of T, changing the flight path from a dive to a 
climb slightly reduces the possibility of oscillatory in- 
stability. 

'The resulting effect of changing the flight path is 
that an airplane will have much less spiral stability and 
may have slightly more oscillatory stability in a climb 
than it will .have in. a dive. This result indicates that, 
for satisfactory stability throughout the flight range, 
an airplane should have more lateral stability than is in- 
dicated by results for level flight. 

As stated before, the computations do not include the 
direct effects' of the propeller nor of the slipstream. 
The large change in power required between a glide and a 
climbing condition is known to have a large effect on the 
stability derivatives of some airplanes. 

Effect of changes in normal acc e leration on lateral 
stability .- Normal acceleration has no effect on the spi- 
ral stability boundary, because it enters into the E 
term of the equations of motion as a multiplying factor. 
(See reference 8.) It increases the magnitude of the 
value of E with respect to the other terms, however, and 
the spiral motion will therefore be affected. 
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Increasing the normal acceleration increases the 
value of Sf/S w required and decreases the maximum value 

of T allowable for oscillatory stability, (See fig. 4.) 
These results indicate that an , airplane which is normally 
static may he unstable when subjected to normal accelera- 
tions. The computations were "based on the assumption that 
the airplane is in straight level flight and is subjected 
to a continuous normal acceleration, This flight condi- 
tion is a physical impossibility. It is believed, however, 
that these flight conditions may bo approached for short 
intervals during pull-ups, in banked turns, and in gusts. 

It was thought that possibly, in practice, the period 
of the unstable motion .of the airplane might be too long to 
be of any particular importance because of the short time 
that the acceleration would exist. Accordingly, computa- 
tions were made for the motions in roll 0, in sideslip 
3, and in heading \|/ with unit rolling moment and with 
unit yawing moment, -each applied independently with normal 
accelerations of lg and 6g. The motions were obtained by 
solving the equations of motion. The assumed characteris- 
tics of the airplane are the same as wore used in obtain- 
ing the stability boundaries. A value for S^/S^ of 0*05 
.and a value for T of 6° was chosen for the computations 
in order that the oscillatory notion would be stable for 
lg acceleration and unstable for 6g acceleration. (See 
fig. 4(c).) 

The resulting notions of the airplane are given in 
figures 5 to 7. In figures 5 to 7 the notions are given 
as- time in seconds to obtain the angle of bank 0, the 
•angle of sideslip P, and the angular heading \|/ in 
radians far unit applied rolling moment or unit appliod 
yaw.ing moment . 

The. motions for 6g acceleration were computed for a 
value of 0^ of 1*4. The motions were also computed for 
lg acceleration at values' of of 1.4 and 0.233 for 

comparisons. .at the same attitude and at the same airspeed. 
In order to obtain an acceleration of 6g, it is necessary 
for this assumed airplane to travel at a speed of 328.7 
foot per second as compared with 134.2 feet por second for 
lg at a value of of 1.4. 

Figures 5 and 6 show largo oscillations in roll and 
sideslip within 3 seconds for the condition with 6'g ac- 
celeration although, in this same time interval, there is 
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no evidence of oscillation with Ig acceleration for the 
same disturbance. These oscillations are likely to con- 
fuse the pilot as well as to allow the airplane to at- 
tain attitudes in which the aerodynamic characteristics 
and the stalling characteristics are greatly different 
from tho se- in"- the normal flight range* 

General comment s The results of this investiga- 
tion are intended to indicate general effects of changes 
in certain airplane parameters and flight conditions and 
arc not intended to indicate the' stability characteris- 
tics of any -particular airplane. The general results of 
the investigation may "be summed up by pointing out the 
factors that make t he ' at t ainment of lateral stability more 
critical and difficult, and the conditions of flight in 
which lateral instability is. most likely to be encountered; 

1, Small amounts of side area for equal values of 0 n ^ 

2* A climbing attitude for spiral instability and a 
dives for oscillatory instability 

3. Any maneuver in which tho normal acceleration is 
increased 

Other factors, such as the type and degree of stability de- 
sired and tho choice, of aerodynamic characteristics used, 
have a large effect on the interpretation of these results. 

Although the typo and tho degree of stability that af- 
fects the control and the riding qualities of the airplane 
*re outside the scope of this report, some' of the follow- 
ing factors should be mentioned; It is normally consid- 
ered that, f or - sat i sf act ory stability characteristics, the 
oscillatory motion should be highly damped, that is, the 
airplane should have a large amount of oscillatory stabil- 
ity. This stability can be obtained by proportioning the 
effective fin area and tho effective dihedral angle in 
order that their values 5 when plotted on a diagram such as 
figure 1 for a particular airplane , will be in a stable 
area and well away from the zero oscillatory stability 
boundary of the airplane. The distance from the boundary, 
although an indication, is not a quantitative measure of 
the amount of damping. For good riding qualities in rough 
air, the amount of spiral stability should be small; some 
stability is considered desirable, particularly for fly- 
ing conditions of poor visibility, although spiral insta- 
bility is generally considered to be preferable to a 



15 



poorly damped oscillatory motion* The value of C n o 
should he positive, giving weathercock stahility; the value, 
however, is dependent upon the , spiral and oscillatory sta- 
hility requirements* These concepts are further substanti- 
ated hy the experimental data of reference 9 in which are 
reported the results of a free— flight— tunnel investigation 
in which the directional stability and dihedral were. varied 
over a wide range 0 In another free— flight— tunnel investi- 
gation (reference* 10 ) it was found that although lateral 
stahility was easier to obtain with larger values of the 
side— force derivative Oy^, , the lateral control character- 
istic a vere considerably impaired and in certain configura- 
tions increasing Oy^ gave poorer over— all lateral flying 

characteristics* . . 

CONCLUSIONS 



From the analysis, for spiral and oscillatory lateral 
stahility in which assumed data are used. for changes in 
certain airplane parameters, the following conclusions may 
bo drawn : 

1 # . An increase in the .aspect ratio of the wing of an 
airplane requires a smaller fin area but makes little dif- 
ference in the attainment of lateral stabil-ity. 

2. Adding side area makes, the attainment of lateral 
stability considerably less difficult. 

3 # An airplane is more likely to bo spirally unsta- 
ble and slightly loss likely to have oscillatory insta- 
bility in a climbing nttitude than in a gliding attitude. 

4* Normal accelerations are likely to- cause oscilla- 
tory instability and the resulting motions may be confusing 
to t he . pilot . 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee, for Aeronautics. 
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